A Cat-State Benchmark on a Seven Bit Quantum Computer 
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We propose and experimentally realize an algo- 
rithmic benchmark that demonstrates coherent 
control with a sequence of quantum operations 
that first generates and then decodes the cat state 
(| 000 . . .) + I 111 . . ))/y/2 to the standard initial state 
1 000. . .). This is the first high fidelity experimen- 
tal quantum algorithm on the currently largest 
physical quantum register, which has seven quan- 
tum bits (qubits) provided by the nuclei of cro- 
tonic acid. The experiment has the additional 
benefit of verifying a seven coherence in a generic 
system of coupled spins. Our implementation 
combines numerous nuclear magnetic resonance 
(NMR) techniques in one experiment and intro- 
duces practical methods for translating quantum 
networks to control operations. The experimen- 
tal procedure can be used as a reliable and effi- 
cient method for creating a standard pseudo-pure 
state, the first step for implementing traditional 
quantum algorithms in liquid state NMR. The 
benchmark and the techniques can be adapted 
for use on other proposed quantum devices. 

Quantum information processing (QIP) offers signif- 
icant advantages over classical information processing, 
both for efficient algorithms ]30| , p9[ and for secure com- 
munication f33|j3f . As a result it is important to estab- 
lish that sufficient and scalable control of a large num- 
ber of qubits can be achieved in practice. There are 
a rapidly growin g n umber of proposed device technolo- 
gies [|J^g||,0@||0] for QIP, and to compare them 
it is necessary to establish benchmark experiments that 
are independent of the underlying physical system. A 
good benchmark for QIP should demonstrate the ability 
to reliably and coherently control a reasonable number of 
qubits. This requires that elementary operations can be 
implemented with small error regardless of the state of 
the qubits, as sufficiently small error is one of the most 
important prerequisites for robust QIP p8|,p^,18,20|]. The 
cat-state benchmark proposed here is perhaps the sim- 
plest demonstration of control which can be implemented 
for any number of qubits and involves coherence in a non- 
trivial way. 

To explain and realize the cat-state benchmark we use 
the example of NMR based QIP. At least two propos- 
als for quantum devices are based on using nuclear spins 
controlled by radio frequency (RF) fields: The first in- 
volves the use of molecules [ [10|Jl4| ] forming an ensemble of 
quantum registers and the second uses nuclei embedded 
in a semiconductor p5|jL7| |. Of these proposals, the first is 



presently accessible to experimental investigation by the 
use of off-the-shelf equipment for liquid state NMR. By 
means of the technique of preparing pseudo-pure states, 
it is possible to benchmark quantum algorithms involv- 
ing up to about ten qubits to determine how well coher- 
ence is preserved and to measure how reliable the avail- 
able control methods are. There have been numerous 
experiments implementing various quantum algorithms 
on up to five qubits P,|,@|@,P,|ll|,||,|| using NMR. 
The experiment reported here coherently implements a 
quantum algorithm on seven qubits with a verifiable fi- 
delity. It also introduces a reliable method for preparing 
pseudo-pure states and for verifying maximal coherences 
in generic spin systems. 

NMR QIP uses spin - nuclei as qubits. Examples are 
protons and carbon 13 bound in a molecule. QIP re- 
quires the ability to couple different qubits. In molecules 
in a liquid at high magnetic field, scalar couplings can 
be used for this purpose and controlled with refocusing 
methods |]lo| , |l^ . Thus, each molecule can be considered 
as a quantum register consisting of (some of) its spin | 
nuclei. The initial state is prepared by allowing enough 
time for thermal relaxation and readout is performed by 
an ensemble measurement using standard NMR methods. 
We use deviation density matrices |HJ for describing the 
state of the nuclei. To simplify the discussion, we use 
a three qubit example of the cat-state benchmark. The 
thermal equilibrium state of a molecule with one pro- 
ton and two carbon 13 nuclei at high field in a liquid 
is given by cr z H) + .2ba z Cl) + .2ba z ° 2) with high ac- 
curacy, up to an overall scale factor and a multiple of 
the identity. The standard Pauli matrices are used as 
an operator basis, and superscripts on operators refer to 
the particle the operator acts on. The cat-state bench- 
mark for this system begins by eliminating signal from 
the carbons to obtain the initial state a z H K Next a se- 
quence of quantum gates Q is used to achieve the state 

o'y 11 ^ (Jy Cl ^ ax° 2 ^ (Fig.|l|). This state is a sum of several 
coherences ]l3|. In particular, it contains the three co- 
herence |000)(111| + |111)(000|, which is the deviation of 
the operator for the cat state (|000) + 1 1 1 1) ) / a/2- If each 
qubit is rotated by a phase <p around the z-axis, the three 
coherence rotates by 30, while all other components of 
this (or any other) state will rotate by 0, (j) or 2<j). This 
feature can be used to label the three coherence and elim- 
inate all other components of the state, for example by 
using a magnetic field gradient [|l3| . An efficient alterna- 
tive using z-pulses or phase cycling is given below. The 
three coherence can be decoded to the state ' |00)(00| 
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(Fig.|T|). This is then observed after inverting the label- 
ing gradient at three times the original strength. In a 
fully resolved reference spectrum obtained by applying a 
90deg rotation to the proton in the initial state, the pro- 
ton shows four peaks, one for each of the states |00)(00|, 
|01)(01|, |10)(10| and |11)(11| of the carbons. After de- 
coding the cat state, only a single peak should be left in 
the spectrum. The ratio F of the intensity of this peak 
to the intensity of the corresponding peak in the refer- 
ence spectrum is unity if everything works perfectly. F is 
reduced by errors in the preparation and decoding steps. 
Under the assumption that error in the phase labeling 
method is negligible, it can be shown that F is a lower 
bound on the average of the fidelities |2(J with which the 
decoding procedure maps the states |000) =b (1 1 1 1 to the 
states (|0) ± |1))|00). 

The three qubit cat-state benchmark can be general- 
ized to any number n of qubits by repeating the steps of 
the cascade in the networks shown in Figjl| We imple- 
mented the seven qubit version using fully labeled trans- 
crotonic acid (FigJ|). The qubits are given by the spin 
| component of the methyl group, the two protons adja- 
cent to the double bond and the four carbon 13 nuclei. 
A fidelity of .73 ± .02 was achieved. The loss of signal is 
primarily due to spin relaxation, incomplete refocusing of 
couplings and intrinsic defects in using selective pulses. 
The success of the experiment derives from the use of 
the following techniques: 1. An RF imaging method to 
greatly reduce the effects of RF inhomogeneities. 2. A 
gradient based selection method for removing signal from 
the spin | component of the methyl group in an almost 
optimal way. 3. The use of abstract reference frames 
for each nucleus to absorb chemical shift and first order 
off-resonance effects in selective pulses. 4. Precompu- 
tation of coupling effects during pulses. 5. A pulse se- 
quence compiler that optimizes delays between pulses for 
achieving the desired amount of coupling evolution while 
minimizing unwanted couplings. All these techniques are 
scalable in principle. Further details are in the methods 
section. 

The cat-state benchmark has three applications that 
promise to make it useful for NMR and other quantum 
technologies. First, the benchmark demonstrates the 
ability to reach the maximum coherence with little loss 
of signal. Previous experiments have generated coher- 
ences by exploiting symmetry and effective Hamiltonian 
methods. Ver y h igh order coherences can be observed 
in solid state JL2J. A maximal coherence of order seven 
was detected by Weitekamp et al. [[32| in benzene with 
one carbon labeled by exploiting symmetry. The meth- 
ods used in these cases do not yield the amount of signal 
that can be achieved by using methods based on quantum 
networks. 

Second, the output of the benchmark can be used as a 
very reliable pseudo-pure state for quantum algorithms. 
We can write the maximum coherence on n qubits as a 
sum of two operators X = |00 ...)(!!.. .| + |H . . .)(00 . . .| 



and Y = z|ll...)(00...| - z|00 . . . . .|. The decod- 
ing operator converts X to |0 . . .)(0 . . .| and Y to 

dy 1 ^ |0 . . .)(0 . . .|. These states can be used as a pseudo- 
pure input to a quantum algorithm using one less qubit, 
provided the following two problems are addressed: The 
first problem is to ensure that the labeling method can be 
used together with a subsequent algorithm. The second 
problem is to eliminate errors accumulated when decod- 
ing the n-coherence. 

Clearly, the method used to label the n-coherence must 
be reliable. The gradient based method is effective in con- 
junction with a quantum algorithm, as long as the echo 
pulse is applied just before the final observation. Unfor- 
tunately, diffusion introduces loss of signal at the gradi- 
ent strengths required when used with long algorithms. 
Also, gradient methods do not easily generalize to other 
proposed quantum devices-an important issue in bench- 
marking. To label the n-coherence one can instead per- 
form 2n + 1 experiments, where in the /c'th experiment, 
the gradient is replaced by explicit pulses that rotate each 
qubit by a phase fa = 27r/c/(2n + 1). If Ok is the expec- 
tation of the observable measured at the end of the /c'th 
experiment, then the value o = ^ k 0/; , e -* 27r/cn /( 2n + 1 ) i s 
non-zero only for signal originating at the n-coherence. 
This technique can be applied in any system where it 
is possible to apply ^-rotations reliably. If the phase of 
applied pulses is highly controllable (as is the case in sys- 
tems controlled by RF or optical fields), instead of ap- 
plying explicit pulses to accomplish the z-rotations, one 
can change the reference frame for each qubit, which is 
equivalent to changing the phase of all subsequent pulses 
and the observation reference phase by — fa. This is 
essentially a phase cycling method for selecting the n- 
coherence Jl^] . We have used both the gradient based 
and this phase cycling method with identical results in 
the crotonic acid system. 

The problem of decoding error can in principle be 
solved by using the maximum coherence directly as the 
input for a (modified) algorithm. However it is not possi- 
ble to obtain a reference signal for the n-coherence with- 
out first mapping it to an accessible observable, which 
can involve a loss of signal. Another problem is that 
it may be inconvenient to use the n-coherence instead 
of the more familiar standard pseudo-pure state. Our 
experiments show that we can decode the n-coherence 
to the pseudo-pure state with no detectable error in the 
observed spectrum. There can be error signal in unob- 
served operators which one would like to eliminate from 
future observation. This can be done efficiently by per- 
forming multiple experiments, each with a random phase 
of Odeg or 180deg applied to qubits 2,3..., a technique 
which is a special case of the randomized methods of [[19[ . 
The number of experiments that need to be performed 
depends on the desired level of suppression of possible 
error signals. N experiments result in suppression by a 
factor of 0(1/ y/N). 

The final application of the cat-state benchmark is as 
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an experiment to test the ability to coherently control 
a quantum system and demonstrate a fully coherent im- 
plementation of a non-trivial quantum algorithm. This 
is a critical issue for scalable quantum information pro- 
cessing, as scalable robustness requires that each opera- 
tion has a maximum error below some threshold (which 
may depend on the types of errors) [p] JT^j20[ | . The known 
thresholds seem to be dauntingly small. Nevertheless, 
interesting small scale computations may be performable 
with much higher error per gate. Thus the ability to 
implement the cat-state benchmark with high fidelity is 
a good indication of what types of tasks can be accom- 
plished in the system at hand. In addition, the decod- 
ing algorithm of the cat-state benchmark is an instance 
of the type of process required to perform fault-tolerant 
error-correction [p8[ , which is believed to be a necessary 
subroutine in any large scale quantum computation. Our 
experiment involved a total of twelve useful two-qubit op- 
erations, so the fidelity of .73 suggests an error of about 
.023 = .27/12 per coupling gate. If this degree of control 
were available in the context of quantum communication, 
it would be close to the known thresholds []5|. 

The realization of the cat-state benchmark given here 
is in an ensemble setting. Most proposals for quan- 
tum devices involve individual systems with pure initial 
states. In these cases the benchmark can be modified by 
replacing the ensemble measurements by repetition to 
infer Ok with sufficiently high signal to noise. The prepa- 
ration step is replaced by a network that directly maps 
the available initial state to the cat state. Note that any 
evaluation of a quantum device involves substantial rep- 
etition, essentially replacing the ensemble measurement 
by an ensemble in time. 

Methods. We used a Bruker DRX-500 NMR spectrom- 
eter with a triple resonance probe for our experiments. 
(The triple resonance probe is normally used for proton, 
carbon 13 and nitrogen 15; we used only the first two.) 
All the equipment used was standard with no specialized 
modifications. The chemical structure of trans-crotonic 
acid is given in Fig.||. Deuterated chloroform was used as 
the solvent. The chemical shifts (at 298K and 500Mhz) 
and coupling constants were experimentally determined 
to within .1Hz by direct analysis of the proton and car- 
bon spectra (see Fig.||). This data was used to design 
selective pulse shapes and times. Only 90deg and 180deg 
rotations were used in the pulse sequence. The hard 
and selective pulses were analyzed by simulation on sin- 
gle and pairs of nuclei and represented optimally as a 
composition of phase shifts, g z g z couplings and an ideal 
90deg or 180deg pulse. The simulation is efficiently scal- 
able, requiring 7(7 + l)/2 two qubit simulations for the 
seven qubit register. This permits elimination of most 
first order errors due to off-resonance and coupling ef- 
fects without using specialized shapes. The computed 
phase shifts were absorbed into the rotating frame of 
each nucleus, while the computed coupling effects con- 
tributed to the coupling operations or were refocused. 
To implement quantum information processing tasks, we 



began with an ideal quantum network expressed in terms 
of 90deg rotations and 1/(2 J) coupling evolutions. Refo- 
cusing pulses are then inserted, and an optimizing pulse 
sequence compiler is used to determine the best choice 
of delays between pulses to achieve the desired evolu- 
tion. The compiler permitted us to automate many of 
the tasks of translating a quantum network to a pulse 
sequence. Often smaller couplings cannot be perfectly 
refocused without an excessive number of pulses. Error 
due to imperfect refocusing is explicitly estimated by the 
compiler. The final pulse sequence used in our exper- 
iment required 48 pulses with an estimated signal loss 
due to coupling errors of .15. 

To have accurate pulses, we have reduced the effect of 
RF inhomogeneities present in standard configurations 
by selecting signal based on RF power. The method first 
applies a 90deg excitation pulse to the methyl group fol- 
lowed by a sequence of pairs of 180deg rotations at phases 
of ±0fc, where was determined from one of the selec- 
tive pulse shapes we used (modified by an initial sequence 
to compensate for an off-resonance effect) and designed 
to cause a 90deg phase shift in the signal at an RF power 
of about ±2% of the ideal. At other powers, the effect 
is such that a phase cycle involving a change of sign of 
the 90deg phase shift eliminates the signal. A final pulse 
returns the selected signal along the z-axis in preparation 
for the next step. By calibrating the power, we were able 
to retain 25% of the signal compared to an unselected 
spectrum. This sequence also has the property of select- 
ing signal from only the methyl group so that the initial 
state is cf z M \ 

The next step in the experiment required selecting the 
spin 1/2 component of the methyl group's state space. 
This can be accomplished by use of a three step sequence 
involving transfer of polarization to the adjacent carbon 
and terminated by a gradient "crusher" (Fig.||). The 
elimination of signal from the spin | states was verified 
by three experiments involving observation of the signal 
on the adjacent carbon after transfer of the methyl po- 
larization with different delays for coupling. One of the 
resulting spectra is shown in Fig Ji| with the standard ref- 
erence spectrum. We were unable to detect error signal 
above the noise. 

The remaining steps of the experiments consist of the 
generation of the n-coherence, labeling the n-coherence, 
and the decoding operations to obtain the standard 
pseudo-pure state, which was observed on the methyl- 
carbon C\ . We chose C\ for making observations because 
all the couplings are adequately resolved there. The se- 
quence is as described earlier, with judiciously inserted 
refocusing pulses and optimized delays. All the pulse 
phases were computed automatically for the nuclei's in- 
dividual reference frames. Knowledge of the intended 
current state of a nucleus was exploited when that state 
is |0)(0| or |1)(1| to absorb the effects of couplings to that 
nucleus into the reference frame. The compiled pulse 
sequence, pulse shapes and other required information 
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needed for running on a Bruker spectrometer is available 
from the authors. Fig.|| shows the pseudo-pure state sig- 
nal compared to a reference spectrum obtained after se- 
lection of the spin | selection sequence on the methyl 
group. Errors can show up as peaks in positions different 
from the leftmost one. We could not detect such errors 
above the noise. The fidelity is given by the ratio of the 
intensity of the left most peak in the final signal to the 
intensity of a peak in the reference spectrum and was 
computed to be .73 ± .02. 
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FIG. 1. Top: Encoding of the deviation matrix a z II 
into a y a y a x by using a cascade of rotations and J-couplings. 
A three coherence |000)(111| + |111)(000| is contained in 
the output which can be labeled using a magnetic gradi- 
ent or phase cycling. Bottom: Decoding the coherence to 
a pseudo-pure state is accomplished by a similar inverse 
cascade. The vertical text below the network denotes the 
state of the three qubits at that point in the network with 
X = o~ x ,Y — cr y , Z — a z . Both networks generalize by ex- 
tending the cascade to more qubits. 
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FIG. 2. Molecular structure of trans-crotonic acid to- 
gether with a table of the chemical shifts and J-coupling con- 
stants. The chemical shifts are on the diagonal and are given 
with respect to reference frequencies of 500.13 MHz (protons) 
and 125.76 Mhz (carbons) on the 500 Mhz spectrometer we 
used. The T2 were greater than 2sec. 
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FIG. 4. Spectra of the methyl-carbon. Spectrum a) is 
obtained after transfer of the methyl equilibrium polarization 
to the carbon. Both the spin \ and the spin § components are 
present. Only spin | signal contributes to the extreme peak 
groups. Spectrum b) is obtained after transfer of the spin 
\ selected polarization from the methyl group. There is no 
signal in the extreme peak groups detectable above the noise. 
This together with two other, similar spectra (not shown) 
obtained after different delays for coupling demonstrates good 
selection of the spin | component of the methyl group. The 
scales for the two spectra are different, the number of scans 
for a) and b) was 8 and 256 respectively. 
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FIG. 3. Network to select the spin \ subspace of the 
methyl group. Refocusings needed to decouple the other 
nuclei are not shown. The spiral at the end is a gradient 
"crusher" to remove the transversal polarization. 





FIG. 5. Spectra of the pseudo-pure state (red) and the 
input state after transfer of polarization to the methyl-carbon 
(blue). Both spectra were acquired with 256 scans and are 
shown at the same scale. The signature of the pseudo-pure 
state is that only a single (the leftmost) peak remains. No 
detectable signal remains in any other peak position, the 
intensity ratio for the left most peak in the pseudo-pure state 
and the reference is .73 =b .02. . 
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